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Abstract 

Using the scheme of mesoscopic nonequihbrium thermodynamics, we 
construct the one- and two- particle Fokker-Planck equations for a system 
of interacting Brownian particles. By means of these equations we derive 
the corresponding balance equations. We obtain expressions for the heat 
flux and the pressure tensor which enable one to describe the kinetic and 
potential energy interchange of the particles with the heat bath. Through 
the momentum balance we analyze in particular the diffusion regime to 
obtain the collective diffusion coefficient in terms of the hydrodynamic 
and the effective forces acting on the Brownian particles. 

1 Introduction 

Brownian effects are ubiquitous in many examples of soft condensed matter 
physics for which the system can be modeUed as a set of interacting degrees of 
freedom in contact with a heat bath. They play a very important role when one 
infers macroscopic behaviors from the mesoscopic level of description, a route 
utilized frequently in the study of complex systems. Many examples, as colloidal 
suspensions polymers H), micelles 0, etc., share those characteristics. 

The dynamics at the mesoscopic level is governed by a set of Langevin pro- 
cesses or equivalently by the corresponding TV- particle Fokker-Planck equation. 
One may arrive at the formulation of those equations by using different meth- 
ods, but the inherent complexity of the description of the dynamics of many 
interconnected degrees of freedom makes it necessary to adopt simple schemes. 

The method we use to analyze the dynamics of a system of interacting Brow- 
nian particles, is based upon a mesoscopic approach proposed recently to obtain 
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the Fokker-Planck equation for the A'"-particle distribution function It ap- 
phes the scheme of nonequihbrium thermodynamics to derive the kinetic equa- 
tion describing the evolution of the A'^-particle probabihty distribution function 

§0. 

One then considers a system of N Brownian particles diluted in a solvent, 
which acts as a thermal bath. The velocity of the particles are viewed as inter- 
nal thermodynamic variables. Such a description, is based on the existence of 
dissimilar time scales for the Brownian particles and permits to proceed with 
an analysis in the phase space of the macroparticles. The local equilibrium 
hypothesis is introduced at the phase space level and from it one derives the 
entropy balance equation. The resulting entropy production accounting for the 
irreversible processes taking place in the phase space allows us to identify the 
fluxes and forces in a similar manner as in the thermodynamic of the irreversible 
processes [0. 

In the linear regime for the dependence of fluxes on the forces one obtains 
the Fokker-Planck equations ||. This theory has been applied to the case 
of interacting Brownian particles under a temperature gradient, and to the 



thermocapillary migration of Brownian droplets |10 . For the case of Brownian 
motion under stationary flow, the mesoscopic approach permits to obtain a 
Fokker-Planck equation which exhibits violation of the fluctuation-dissipation 
theorem i.e., the diffusion coefficient does not obey the Stokes-Einstein 
law. In addition, this formalism has been applied to obtain kinetic equations 
for polymer dynamics [ p2[ . 

The above approach offers a mesoscopic framework for which simply by 
applying the scheme of nonequilibrium thermodynamics in phase space it is 
possible to derive kinetic equations of the Fokker-Planck or Smoluchowski types 
to analyze the dynamics. Our goal in this paper is to apply that method to 
describe the dynamics of interacting Brownian particles in the approximation 
of pair additive interactions. 

The paper has been organized in the following way. In section 2, we perform 
a contraction of the description of the A^-body problem to two tracer particles, 
where their velocities are the internal degrees of freedom from the thermody- 
namic point of view. Section 3 is devoted to obtain the Fokker-Planck equations 
for the one- and two particles distribution functions. In section 4 we discuss the 
hydrodynamic regime, whereas in section 5 we analyze the diffusion regime 
from the momentum balance equation. Finally, in section 6, we present some 
concluding remarks. 

2 Conservation laws for the probability distri- 
bution functions 



We consider a set of N interacting Brownian particles of mass m suspended 
in a solvent which plays the role of a heat bath at rest. We associate to each 
Brownian particle a phase coordinate vector x; that denotes its position r^, 
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and its velocity in phase space, i.e. — (rj,Ui). Hence, the phase space 
occupied by the N particles is the one spanned by the values of the vector 
r — (xi, X2, ...xat) = (x"). The probability density for the A'^ particles to be 
at point r at time t is denoted by p(^)(r,f) and satisfies the normalization 
condition, 

y pW(r,t)dr = 1. (1) 

The entropy of the Brownian particles is given through the Gibbs entropy 
postulate, 

5 = - fc / f(^) (r, t) In ^'^^l^'^^ dr + S,,. , (2) 

Peg. (r) 

where Seq. is its equilibrium value, k is Boltzniann's constant and Pe^'^ (F) 
the iV-particle equilibrium probability density. In the framework of mesoscopic 
non-equilibrium thermodynamics ^ , the phase coordinates of the set of N 
Brownian particles play the role of internal degrees of freedom. The Gibbs 
equation is formulated as follows, 

SS = -^J fi{r,t)Sp(^HT,t)dT, (3) 

where /i(r, t) is the non-equilibrium A^-particle chemical potential, T is the heat 
bath temperature which we assume constant and 6 stands for the exact total 
differential. The rate of variation of the entropy is then given by 



dt T 



y'A*(r,i)|pW(r,i)dr. (4) 



Quite generally, we may assume that p(^)(r, t) satisfies the continuity equa- 
tion, 



OpiN) ^ 9p(N) ^ QpiN) N 

i—l ij = l i—1 

where 0^ represents the direct interaction potential between particles. The 
energy and momentum dissipative interchange between particles and the solvent 
is taken into account through the fluxes ji^'' — j1^''(x^). 

From the A^-particle continuity equation (5), after integrating over N-s phase 
coordinates, we obtain the reduced continuity equations. 



(6) 
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with PW = pW(x^^) 

AT] r{N-s) 

F(^)(xi,X2,...x„^)= ^^J^^J pWdx,+i...dxjv, (7) 
the reduced distribution function and 

Jif = / Ji^^rfx,+i...dx^, (8) 

the reduced fluxes in phase space. 

In particular, we are interested in the case of the one- and two-particle distri- 
bution functions, hence from the above equation we can obtain, the continuity 
equations, 

dt 5ri ""7 dr, 5ui ^ui "^"^ 

and 

1=1 *)j=i *=i 

These cqiiations present the usiial drift terms and conservative interaction con- 
tributions on the left hand side. In eq. (9) the third term expresses the mean 
force over a Brownian particle due to the presence of other particles. On the 
other hand, eq. (10) includes the effect of the direct forces felt by a pair of 
particles. The right hand side of each equation accounts for the dissipative in- 
teraction of the particles with the solvent through the corresponding currents. 
At this stage, we have provided the continuity equations for the reduced proba- 
bility densities P^^^ and P^^' , where the fluxes J^^' appear as unknown functions. 
One of the features of the mesoscopic nonequilibrium thermodynamics approach 
presented here, is to find out the expressions of J^^^ and Jut , which will be 
obtained from the entropy production in the next section. 



3 Fokker-Planck equations for the one- and two- 
particle distribution functions 

The simplest model of N-interacting particles can be implemented by only 
considering interactions between pairs. This approximation implies a substantial 
simplification which in terms of our mesoscopic description concerns the analysis 
of the dynamics through the one and two-particle distribution functions. In 
this section we will derive the kinetic equations of the Fokker-Planck type, 
accounting for the evolution of those distribution functions. 
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To proceed, we factorize the N-particle distribution function P(^^(x^,f) as 
follows 

pW(x^,i) = P(i)(xi,i)p(i)(x2,t)...p(i)(xjv,t)5(^) (x^,t) , (11) 

where P^^^ (xi,t) represents the reduced distribution function for s = 1; more- 
over, g*^^) (x^ ,tj is the A/^- particle dynamic correlation function. In a similar 
manner, as for the equilibrium situation [13-15], we propose a factorization for 
the correlation function g'^^^ {^^^t) = 9^^\^i,^2, •••,xjv,t) as follows, 



{'^N-2,XN-l,^N,t) ...5g^ ^ (Xi, ...,Xjv,t) . (12) 
Here, 5g^^^ (xi,X2,X3,f) is defined through the relation 



(Xi , X2, X3, i) = 5(2) , X2, i) g^^^ (Xi , X3, t) 5(2) (X2, X3, t) 5g^^^ (Xi , X2, X3 , t) , 

(13) 

such that when 6g^^^ (xi,X2,X3,t) = 1, it reduces to the so-called Kirkwood's 

superposition approximation, and g^**) is defined similarly, through Sg^''^\ In the 
above factorization for c/(^\ we have (1/2) iV (iV - 1) pairs, {1/3\)N {N -1) 
{N — 2) triplets, etc. so that 



InPW = ^ln(pW(xi,i))+lnffW(x^,i) 



= 7VlnP(i) (xi,t) + ^{N- l)ln5(2) (xi,X2,t) + 

^ (TV - 1) (iV - 2)ln^5(3) (xi,x2,X3,i) + ... (14) 
Similary, for the equilibrium case [13-15], we have 

InPif ) (x^) = iVlnP^W (xi) + y (iV - 1) \ngf^ (n, r2) + .... (15) 

Substitution of these two last equations in (2), yields the expression 
of the Gibbs entropy postulate. 



PW (xi,f) In-^ 



P« (xi,t) 



dxi 



. + Seq. (16) 
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It is worth mentioning that a similar expression, for the non-equihbrium en- 
tropy, has been obtained for a kinetic description of a dense gas Q and its 
equihbrium counterpart has also been used for a calculation of the thermody- 
namic entropy In addition, this statistical entropy expression has been 
recently reformulated in terms of the potentials of mean force between parti- 
cles and their graphical representation has been presented for the excess 
entropy [ pl[ . 

Similarly to the A^-particle case, we formulate the Gibbs equation as follows. 



SS = -i^|^a)(xi,i)5pWdxi 

where + fi^'^^ represents the non-equilibrium chemical potential per unit 
mass of the suspended particles. Their general forms are given by, 

^(i)(xi,t) = — lnP(i)(xi,t) + ¥^\^it) (18) 

and 

kT 

/i(^)(xi.X2, t) = — ln,g(2)(xi X2, + *(^)(xi X2, i), (19) 
m 

where ^''•^^ and are unknown potential functions. 

At equilibrium, the distribution functions must satisfy the next expressions, 

n^,^.^M = exp^(M^^^^-|), (20) 

where we have defined the distribution function in a reference frame moving 
with the mean velocity of Brownian particles and 

5g) (ri,r2) = exp^ (/if;^ (ri,r2) - cj,,, (ri,r2)) . (21) 

Here /z^^^ and iJ.'g^^ (ri,r2) correspond to the ideal and excess parts of the 
chemical potential of the Brownian particles at equilibrium with the solvent. On 
the other hand, according to the above definition of the pair correlation function 
we can identify the effective pair interaction potential, namely, 4>l2^ (ri,r2) = 
(j)i2 (ri, r2) — /x^'^ (ri, r2) , where —^'j^'^ stands for the indirect part of the work, 
corresponding to the free energy change, to transport the Brownian particles at 
a certain separation, accounting for the static interaction of the particles with 
the solvent. For the particular case when the effective interaction (ri,r2) 
has spherical symmetry, geq only depends on the distance separation modulus 
of the two particles ri2 =| ri2 | = | r2 - ri |, i.e. gig* (| r2 - ri |) 
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If wc now substitute (20) and (21) in (18) and (19), such that in this hmit 
case /x^^) = /i^^^ and fi^^'' = fi'^g^^ , and using the definition of the total chemical 
potential at equilibrium = /^^^ + fi'^^ , one obtains the expressions 



= I, (22) 



and 



= (23) 
With this identification, we can rewrite expressions (18) and (19) as follows. 



and 



M^'^ = — (24) 



m''^ = ^1M^1+Mb:- (25) 



m 



(2) 



After substitution of these expressions in the entropy functional (17), using 
the definition of mass density = m J f^^^dui for the Brownian particles and 
considering that the heat bath density is constant, we can identiiy through 
the Gibbs entropy postulate (16) the equilibrium entropy, 

Mb* 

SSe, = -^Sp^. (26) 

This means that eq. (17) with (24) and (25) is consistent with the Gibbs entropy 
postulate (16) formulated at the pair correlation level. 

We now proceed to analyze the evolution in time of Gibbs local entropy 
S{ri,t). To reach this purpose, from eq. (17) we obtain the rate of change of 
entropy per unit volume. 



|5(r„t) = -miV J ^|pa)rfu,-m^^(Ar-l) J ^|p(2)rfu,rfu,rfr,. 

(27) 

Working out the first term in this expression, after substituting eq. (9) we have, 
. d . d 

m 



J ^§iP'''du, = -^^-ji^\r„t)+a^'\v„t). (28) 



Here 



jW(ri,t) = iVfc|uiP«(^/zW-l)dui 
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2T J \ dr^^ J ri2 7o ^Ui 

(29) 

is the first contribution to the entropy flux, with P^'^^da = P^"^^ (ri — [1 — a] 
rij, Ui, ri + arij, U2,t) da, as results from assuming that P^^-* is a slow varying 
function of ri and it admits a Taylor's expansion around the distance separation 
between particles Vij The first contribution to the entropy production is 
then given by, 

a(i)(ri,t) = -Nk J J« . ^In (^^^ dui. (30) 

We now analyze the two-particle term in eq. (27) using eq. (10). After some 
algebra we obtain the next identity 



where 



2T J 2ri2 dri2 Jq 

(32) 



2r 



is the two-particle entropy flux and 



= ^ 






d 
dui 


In j 




N{N - 
2 




f 


d 

8X12 


In j 


\9if 



duidu2dr2 
duidu2dr2. (33) 

y9eq' J 

is the corresponding two-particle entropy production. 

A previous analysis on simultaneous motion of A^-Brownian particles leads 
to the next global entropy production 



d^'\ (34) 
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Inserting in this expression the distribution functions InP*^^^ and In Peq \ given 
through eqs. (14) and (15) we obtain an alternative expression for the global 
entropy production, i.e., 



/ 8 f p(iM 



N{N-1) 



A; V / Jif • ^ In ^ du2dri2dui + ... (35) 

\9eg.J 



Where we have taken into account the definition of the reduced fluxes, eq. (8). 

We notice that the local entropy production a{ri,t) = a*^^^ (ri , t) + a^^^ (ri , t) 
can be identified in the integrand of eq. (35). This means that consistently we 
can use the integrand of the eq. (34) as the local entropy production, and after- 
wards the hypothesis that the fluxes are coupled to the thermodynamic 

forces ^ In ^^rwy^ ) giving rise to the next phenomenological equations, 

j^ij=l \^eq / 

Defining the friction coefficients as 

= ^^^i^j (-37-) 

we can give an alternative and more useful expression for the fluxes , namely, 

^ ' kTeq aP(^) 



m duj 



(38) 

Particularly, for the one and two-particle description the corresponding fluxes 
are, 

kTea apw\ 



and 

m dui J 



i,.7 = l ^ 



(39) 



(40) 



Inserting now these expressions into the continuity equations (9) and (10) 
one obtains the desired Fokker-Planck equations. 
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and 



— P(i) + ui • ^P^^^ - m W ^fi^ • — du2(ir2 

at OTi 7 OTi aui 



-p(2)+^u,.— rn-i ^ 

^-^ ovi ^-^ ovi aui 

The analysis of these equations, deserves some comments. The third term 
on the left hand side of eq. (41) represents the local equilibrium mean effective 
interaction potential applied at the point ri, where we have included the in- 
stantaneous interaction of a Brownian particle placed at ri, with the solvent. If 
we do not take into account this mean field term, coming from the interactions 
of a cloud of Brownian particles around ri , including the local structure of the 
solvent, we recover the usual one-particle Fokker-Planck equation ||2^ . 

With respect to eq. (42), the effective interactions —m~^ j=i ■ ^g^.' 
are applied in a non-local manner over the two tracer particles placed at ri and 
r2, respectively. The friction tensor is taken by pairs now with respect to the two 
selected particles. This Fokker-Planck equation coincides with the one obtained 
first by Mazo ||2^ and more recently by Piasccki et al This last one was 
derived for a system of two-Brownian hard spheres immersed in a hard-sphere 
solvent. 

When the velocities of the Brownian particles thermalize the system reaches 
the state of local equilibrium. For this case, the dynamical description of the 
Brownian particles would be represented by the evolution equations for the one- 
and two- particle local equilibrium distribution functions. Hence, the right hand 

sides of the continuity equations ( 9) and (10) should be replaced by — j and 

2 9J<^' 

dr' ' respectively. Following a similar analysis as the one performed 

in this section, we would obtain the expressions for the fluxes J^i^ and Jr?''. 
The final result would correspond to a couple of equations of Smoluchowski 
type. Such approach has allowed the calculation of the long time self-diffusion 
coefficient of hard sphere suspensions and, the corresponding comparison with 
Brownian dynamics simulations | |26[ |. 

The one- and two-particles Fokker-Planck equations describe the dynamics of 
the system at mesoscopic level, under the assumption of pair interactions. These 
expressions can be used to analyze the hydrodynamic regime of the Brownian 
gas. This will be the goal of the next section. 
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4 Hydrodynamic equations 



Our purpose in this section is to derive the complete set of hydrodynamic 
equations describing macroscopicaUy the dynamics of the gas of Brownian parti- 
cles. For simplicity in our notation, we will omit the dependences of the different 
quantities in (ri,i). The conservation law for the mass of the particles follows 
from the definition of the density 



F(i)dui. (43) 



We take the temporal derivative of this expresion and use the continuity equation 
(9). After integrating by parts the resulting equation and using the fact that 
the current 3^} decays very rapidly when increasing the velocity, one obtains 
the conservation law 

We can proceed in a similar wayQ with the definition of the momentum 
density of the Brownian particles 

p^WB^mJuiP'-'Uui. (45) 

By taking the time derivative of this equation and using the one-particle Fokker- 
Planck equation (41), we obtain a preliminar balance equation for the momen- 
tum density. 



Here 



- J P('^^dr2du2dui - /3nP,VB. (46) 

y|=mypW(ui-VB)(ui-VB)dui, (47) 

is the kinetic part of the pressure tensor for the Brownian particles, with 
the friction coeficient, such that PuP^'Vb takes into account the momentum 
exchange between the Brownian particles and the bath. 

Now, we admit that the distribution function P^^^ varies slowly with the 
space coordinates as it should be for a concentrated suspension. Expanding the 
probability density P^^^^ around the distance ri2 = r2 — ri between particles 
[ETI, we may approximate 



/P(^)^dr2 
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~ -T. — TT / ri2ri2^2 / P'^^^ (ri - [1 - a] ri2, ri + ari2, Ui U2,t) dadri2, 

ori 2 J ri2 Jq 

(48) 

with 0]^2 ('"12) — §777 ■ This relation enables us to identify the potential compo- 
nent of the pressure tensor, namely, 



^ ' ; ^12 ('^12) 
V B = -TT / ri2ri2-t^-^ X 



27 ri2 

y P"^^^ (ri - [1 - Q;]ri2,ri + ari2,Ui,U2,t) da^ dri2duidu2. 

(49) 

By means of the above identity, the momentum balance (46) can be rewritten 



in the usual manner 123, 28 



with 

^B = 'v'h + 'VB, (51) 

the total pressure tensor of the Brownian particles and ^ = ^ + v_b • ^ the 
hydrodynamic derivative. This expression for the pressure tensor is consistent 
with the one obtained by Felderhof j2^ (eqs. 3.11 and 3.14 of such reference), 
who used the N-particle Fokker-Planck equation for its derivation. 

In order to complete our scheme, we need to obtain the balance equation for 
the density of internal energy of the "gas" of Brownian particles . This 

quantity splits up into kinetic and potential contributions 

where 

= f j P''Hn,-.sfdu, (53) 

and 



4> m 
p u — — 

^ B B O 



j P^^Ui2duidu2dr2, (54) 



are the densities of kinetic and potential energies, respectively. Taking the 
temporal derivative of eq. (53) and using the Fokker-Planck eq. (41), one 
obtains, 
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d k dm 



-•- / (ui - VB)ViPW(iui + y / pWui- 



9ri 



m oui 



d 
Or I 



Defining now 



J (ui - vb)' (ui - vb) P^^^rfui + ^ / ('^i - ^s)' -P^'^'^uivs 



-m (^^) / P^'^ (ui - vs)^ dui. 
= f / (ui - vs)" (Ui - Vs) P^'^dui 



(55) 



the kinetic part of tlie lieat flux, using tlie deflnition of the kinetic energy density 
(53) and the next identity, 

P^^^ui ■ (ui - vs)^ dvLi = 



= m 'ui (ui — vs) : 5 "Ui 



dri 



m 



j (Ui - Vb)(ui - VB)(iUi 



dvi ' 



we obtain the balance equation for the kinetic part 



d k 

at'-''- 



dri dri 



(56) 



The derivation of the above balance equation deserves some comments. First 
of all we have assumed that the friction coefficient is independent of the veloc- 
ity. This means that we have not taken into account the non-linear friction 
case. The velocity dependence of the friction coefficient (3^^ (ui) concerns with 
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the case in which the Brownian particles are able to take up external energy, 
which can be stored in an internal energy depot This phenomenon is the 
so-called active motion and is of interest in the dynamical analysis of driven 
physico-chemical and biological systems If we want to take into ac- 

count this fact, the above balance equation (56) would include the additional 
term — m / p(-^^/3^^ (ui — VB)-VB(iui. Moreover, the temperature T corresponds 
to a non-equilibrium distribution of the Brownian particles, which eventually 
thermalizes to the usual Maxwell velocity distribution function Pig^with the 
equilibrium temperature T^q. In this manner, the last term of eq. (56) accounts 
for the kinetic energy interchange between the Brownian particles and the heat 
bath. 

Using the kinetic energy density definition (53) together with the energy 
equipartition theorem, namely, 

= ^^^T^eg, (57) 

the last term in the equation (56) can be rewritten as follows. 



- - (^^) / P''' K - du, = -2pJ^^ - <^ ) . (58) 

With this expression, we recover the usual kinetic energy balance equation 
for a dilute suspension of Brownian particles H , 



The potential energy balance equation, follows from its definition eq.(54). 
By using the two-particle Fokker-Planck eq.(42) and Newton's third law ^q^. = 
— we obtain 



- • + PsuI-b) + f / P(^) (u, - u,) . ^du,du.dr„ 

(60) 



in which 



(n, i) ^mjiui- vs) ^P^^^duidu2dr2, (61) 

is the current of potential energy density transport and ^ J P^^^ (ui — U2) • 
^S^duidu2dr2 is the production term. After applying the approximate expres- 



sion (48) one has 



O <t> d /,(!) ,(2) * 
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./P(^)|li.u,du,«r.-^^|^, (62) 



where 

<\>\2 ('"12) 



■r\2 



Ui + U2 - 2vb) X 



X 



/ mP(^^ (ri - [1 - a] ri2,U2,ri +ari2,ui;t)daduidu2dri2 (63) 
Jo 



is the potential contribution to the heat transport. 

Regarding the term J P*-^-*^^ • Ui(iui(iu2(ir2, when we neglect the veloc- 
ity dependence of the pair dynamic correlation function, i.e., P^^-* (xi_X2, t) = 
P(i) (xi, i) P(i) (x2, t) (n r2,i), we have, 



/ P(2)^.uiduidu2dr2 

J OTi 



VB • / ^' g^^^ (ri r2.t) ^dr2 
m J m ori 

= -'-^^.,.Fir,,t). (64) 
m 

Taking into account that the potential contributions to the pressure tensor 
(eq. 49) and the second part for the heat flux (eq. 63) estimate the short 
range interactions by means of the approximation described in eq. (48), the 

force F(ri,f) = — / g^^'' (ri,i'2.^) ^^f-dr2 accounts for the long range 

interactions. With these considerations, the potential energy balance takes the 
conventional macroscopic form , 

dlP.^s - -d^, ■ W + + P^^s-b) -'Pb--^---b-F. (65) 
The balance of internal energy follows just by adding eqs. (59) and (65), i.e, 
qIPb^b + • (J, + Pb^s-b) = -F+P^: — 

-2pB/3,, (<-<J, (66) 

where 

Jf =V- + JS+JS (67) 
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is the total heat flux, due to the energy interchange of Brownian particles with 
the solvent. The term — — vr • F accounts for the rate at which work is done 

m ^ 

by a Brownian particle over the solvent due to the static pair effective long 
range interactions. This term would be useful to account for the active motion 
mentioned above pO[ |. 

In this manner, with the help of two-particle Fokker-Planck eqs. (41) and 
(42), we have obtained the corresponding balance equations for the suspended 
particles. The expressions for the pressure tensor (eqs. (47) and (49)), and the 
heat flux (eqs. (55), (61) and (63)) are of particular interest. Using a suitable in- 
teraction potential and a pair correlation function for the interacting Brownian 
particles, and comparing with Fourier and Newton laws, approximate relation- 
ships for the heat conductivity and viscosity can be obtained. Moreover, with 
the help of the momentum balance eq. (50), we can analyze the diffusion regime. 
This point will be the treated in the next section. 



5 Diffusion regime 

As follows from the balance equation for the momentum of the Brownian 
particles, the mobility bn = , introduces a characteristic time scale defining 
the inertial regime in the dynamics of the Brownian particles. For times such 
that t » the particle enters the diffusive regime in which the inertial term 
on the left hand side of eq.(50) is negligibly small, hence this equation yields, 

-^•"^B^/SnP^VB. (68) 

For the particular case where no external flow is imposed to the system, the 

pressure tensor V b reduces to the osmotic pressure V b — Pb ^ j with U 
being the unit tensor and 

the pressure that the set of Brownian particles exerts on the host fluid, recover- 
ing the usual virial equation of state. This result agrees with the one obtained 
by Brady using purely hydrodynamic arguments | |3^ . 

Using the definition of the mass current J = p^v^ and the expression of the 
pressure (69), we obtain from eq. (68) the Fick's law 

J = -D,^, (70) 
from which we can identify the collective diffusion coefficient. 

kT ' 



hi. (71) 



For long times, this quantity accounts for the relaxation of the spatial distri- 
bution of Brownian particles. In the redistribution of particles, we have included 
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the influence of the solvent through the additional contribution to the potential 
interaction per particle. For the particular case in which the effective interac- 
tions related to (f)i2 decay faster than r^^ , where ri2 is the relative separation 
of two Brownian particles, the diffusion coefflcient becomes 



Dc^ — 

TO 



l-47r^ 







bii, (72) 



which agrees with the corresponding expression given by Russel [^5| . 

Alternatively, by using the virial equation of state for the Brownian particles 
(69), we obtain the usual thermodynamic relationship [|3|, Q, between the 
macroscopic collective diffusion tensor and the osmotic compressibility, 

The ensuing expression relating collective diffusion, mobility and compressibil- 
ity, deserves some additional comments. First of all, it relates the mass transfer 
ratio with the isothermal compressibility as the driven force and the hydrody- 
namic mobility as the response function. The mobility appears in our scheme as 
a phenomenological coefficient, which may in general depend on timcp^. The 
explicit expression of this coefficient must be borrowed from hydrodynamics. 

At the two-particle level of our description, the appropriate approach is the 
one given by Batchelor Moreover, the pressure tensor (51) arising from 

the momentum balance accounts for the pressure exerted over the solvent, thus 
the pressure difference between two separate points drives the collective motion 
of Brownian particles from different sides of the solvent. This fact has been 
used to interpret light scattering experiments in micellar emulsions Q , Q] and 
lyophilic silica particles in nonpolar solvents where eq. (71) and the virial 
form for the pressure were used as a working hypothesis. It has served to test 
the application of liquid theories to suspensions and to analyze its validity for 
arbitrary concentrations and for pair continuous interaction potentials ]35[| [^6| . 

The above relation must be complemented with the functional dependence 
of bii ($) and pb {^) , on $ = ^na^^ the volume fraction of the Brownian 
particles. The expression for the mobility bn ($) can again be obtained from 
hydrodynamics. For a hard sphere model with stick boundary conditions in the 
dilute regime, the result is bn (<i>) — (Gnrja) ^ (1 — 6.55$) with 77 the vis- 
cosity of the fluid phase and a the radius of the Brownian particle. An analysis 
that makes use of a combination of stick and slip boundary conditions helps 
to fit experimental data|^. Moreover, for strongly charged spheres at low salin- 
ity based on an effective macroion fluid theoretical model, the mobility is well 
represented by the parametric form bn ($) — {6Tn]a)~^ (1 — p^") Using a 
model of effective hard spheres with $ -dependent diameter, the values p ~ 1.8 
and a = |, can be explained. 

The effect of interactions in the Brownian motion, was discussed in a phe- 
nomenological manner by Van den Broeck et al |Q through the generalized 
Einstein expression (73). They use the Batchelor's expression for the mobility 
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and the virial form for the osmotic pressure (69), such that for the case when 
the interaction potential contains in addition to a hard core part a purely at- 
tractive interaction the mobility bn (<i>) increases and both, the compressibility 

(ip^) ^^^"^ collective diffusion Dc decrease. The opposite effect occurs 
when a purely repulsive interaction appears along with the hard repulsive part. 
These results have been corroborated by experiments ||^ . 

We would like to notice on the fact that although our study has been per- 
formed for interacting Brownian particles the obtained collective diffusion coef- 
ficient (73) could be used to study sedimentation processes j36|. Moreover, to 
interpret recent experiments of sedimenting particles , the functional form of 
the osmotic pressure pB (*&) for hard sphere Brownian particles is particularly 
useful. Our study focused on continuous pair interactions, we then confirm the 
validity of the relation between collective diffusion, mobility and compressibility 
(eq. 73) for any pair effective interaction. 



6 Concluding Remarks 

In this paper we have analyzed the dynamics of a set of interacting Brownian 
particles by using the method of nonequilibrium thermodynamics. Our approach 
starts from a previous treatment of simultaneous Brownian motion of colloidal 
particles, for which the Fokker-Planck equation, describing the evolution of the 
TV-particlc distribution function, was derived in the framework of mesoscopic 
nonequilibrium thermodynamics. The evolution equations for the s-particle 
distribution functions arc obtained from the continuity equation in phase space. 
The currents ocurring in those equations follow from mesoscopic nonequilibrium 
thermodynamics. The first two equations are the one- and two-particle Fokker- 
Planck equations. 

The corresponding Fokker-Planck- like equations are analogous to the ones 
obtained by Mazo j2J] and Piasecki et al. j2^, for the case of hard-sphere 
macroparticles in a hard-sphere fluid with a multiple time scale which permits 
the derivation of microscopic expressions for the friction tensors between the 
macroparticles. The account of the discrete nature of the solvent is very useful 
to describe, for example, the dynamics of nanocoUoids ]44| . 

From the kinetic equations we have discussed the hydrodynamic regime. The 
balance equations for the conserved quantities are derived from the equations 
relating those quantities with the moments of the distributions. The momentum 
balance permits the analysis of the diffusion regime emerging at long times. In 
this particular level of description, we obtain a generalized Einstein relation for 
the collective diffusion, mobility and compressibility of the Brownian particles, 
where the excess chemical potential and the mobility stand for the fit parame- 
ters. 
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